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Particle decay during inflation is studied by implementing a dynamical renormalization group 
resummation combined with a small A expansion. A measures the deviation from the scale invariant 
power spectrum and regulates the infrared. In slow roll inflation, A is a simple function of the slow 
roll parameters ev, rjv- We find that quantum fluctuations can self-decay as a consequence of the 
inflationary expansion through processes which are forbidden in Minkowski space-time. We compute 
the self-decay of the inflaton quantum fluctuations during slow roll inflation. For wavelengths 
O , deep inside the Hubble radius the decay is enhanced by the emission of ultrasoft collinear quanta, 

• i.e. bremsstrahlung radiation of superhorizon quanta which becomes the leading decay channel for 

|_J ' physical wavelengths H <C k v h{jj) <C H/(t]v — ev)- The decay of short wavelength fluctuations 

hastens as the physical wave vector approaches the horizon. Superhorizon fluctuations decay with 
a power law rf in conformal time where in terms of the amplitude of curvature perturbations A^, 
the scalar spectral index n B , the tensor to scalar ratio r and slow roll parameters : 

<N ' 
O i 

, The behavior of the growing mode rp v Ev jr\ features an anomalous scaling dimension F. We 

■ discuss the implications of these results for scalar and tensor perturbations as well as for non- 

f"-) ' gaussianities in the power spectrum. The recent WMAP data suggests F > 3.6 x 10 — ° 

o 

^3 

6 

+3 . A period of accelerated expansion in the early universe, namely inflation, is nowadays part of standard cosmology 
since explains the homogeneity, isotropy and flatness of the observed Universe At the same time, inflation pro- 

vides a mechanism for generating metric fluctuations which seed large scale structure: during inflation physical scales 
grow faster than the Hubble radius but slower than it during both radiation or matter domination eras, therefore 
physical wavelengths cross the horizon (Hubble radius) twice. Quantum fluctuations generated during inflation with 
wavelengths smaller than the Hubble radius become classical and are amplified upon first crossing the horizon. As they 
re-enter the horizon during the decelerated stage these fluctuations provide the seed for matter and radiation inho- 
mogeneities which generate structure upon gravitational collapse 0, H> 0> El • Most of the inflationary models 
predict fairly generic features: a gaussian, nearly scale invariant spectrum of adiabatic scalar and tensor primordial 
perturbations (gravitational waves). These gene ric predictions are in spectacular agreement with Cosmic Microwave 
Background (CMB) observations. Gaussian |13| and adiabatic nearly scale invariant primordial fluctuations ^3 pro- 
vide an excellent fit to the WMAP data as well as to a variety of large scale structure observations. Perhaps the most 
striking confirmation of inflation as the mechanism for generating superhorizon ('acausal') perturbations is the anti- 
correlation peak in the temperature-polarization (TE) angular power spectra at I ~ 150 corresponding to superhorizon 
scales [TEL Hi! . The anticorrelation between the E-mode (parity even) polarization fluctuation and the temperature 
fluctuation is a distinctive feature of superhorizon adiabatic fluctuations |17| : the (peculiar) velocity gradient gen- 
erates a quadrupole temperature anisotropy field around electrons which in turn generates an E-polarization mode. 
By continuity, the gradient of the peculiar velocity field is related to the time derivative of the density (temperature) 
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fluctuations, hence the peculiar velocity and the initial (adiabatic) contribution to the (acoustic) oscillations of the 
photon baryon fluid are out of phase by 7r/2 [17|. Thus, the product of these two terms gives an anticorrelation peak 
at kc s rjdec = 3 7r/4 which corresponds to super horizon wavelengths since the size of the horizon is y3 larger than the 
size of the sound horizon. The WMAP (TE) data [TElllr3 | clearly displays the anticorrelation peak at I ~ 150 providing 
perhaps one of the most striking confirmations of adiabatic superhorizon fluctuations as predicted by inflation. While 
the robust predictions of a generic inflationary model seem to provide an excellent fit to the WMAP data, different 
models predict slight differences. Therefore, theoretical differences between different models as well as potential exper- 
imental deviations from the most generic features are the focus of intense study. With the ever increasing precision of 
CMB observations it is conceivable that forthcoming observations will allow a narrower determination of inflationary 
models. Relevant discriminants between models are: non-gaussianity, a running spectral index either for scalar and/or 
tensor perturbations, an isocurvature component of primordial fluctuations, etc. Already WMAP re por ts a hint of 
running spectral index of scalar perturbations from the blue on large scales to the red on small scales |lfij . Quantum 
effects associated with interactions can potentially lead to non-gaussian correlations^^ - !!!!- Therefore the detection 
of a running index (as hinted in the WMAP data) or small non-gaussianities in the temperature correlations imply 
potentially interesting quantum phenomena during the inflationary stage that was imprinted on superhorizon scales. 

The inflaton is usually studied as a homogeneous classical scalar field0, 0, @|. However, important aspects of the 
dynamics require a full quantum treatment for their consistent description. The quantum dynamics of the inflaton is 
systematically treated within a non-perturbative framework and some consequences on the CMB anisotropy spectrum 
were analyzed in ref . |24j . 

In this article we study quantum phenomena during inflation which contribute to relevant observables in the CMB 
anisotropics and polarization. In particular, we focus on inflaton decay during inflation as a potential source of 
quantum phenomena contributing to deviations from scale invariance in the primordial power spectrum and/or to 
non-gaussian features. If the inflaton couples to other particles, then its quantum fluctuations which seed scalar 
density perturbations also couple to these other fields. Consequently, the decay of the amplitude of the quantum 
fluctuations of the inflaton may lead to a modification of the power spectrum of density perturbations. The same 
coupling that is responsible for the decay of the inflaton quantum fluctuations can be also the source of non-gaussian 
correlations. 

Particle decay is a distinct feature of interacting quantum field theories and is necessarily an important part of 
the inflationary paradigm: the decay of the inflaton into lighter particles after inflation may yield to the radiation 
dominated stage. Recently, inflaton decay during a post-inflationary stage has been considered as a possible source of 
metric perturbations arising from an inhomogeneity in the inflaton coupling 25] . Inflaton decay has also been studied 
as a dissipative mechanism in the dynamics of the inflaton 26] , however these studies only apply when the expansion 
rate is much smaller than the typical mass scales. 

In a previous article |27| we introduced and implemented a systematic program to study the relaxational dynamics 
and particle decay in the case of a rapidly expanding inflationary stage. Whereby rapid expansion refers to the Hubble 
parameter during inflation being much larger than the mass of the particles. In the case of the inflaton, this is the 
situation of relevance for slow-roll inflation and a necessary (although not sufficient) condition for an almost scale 
invariant power spectrum of scalar fluctuations 0, 0, IE H3 • In ref. [2£| inflaton decay was studied in some particular 
cases for which a solution of the equations of motion was available. The method of ref. [2^] was recently applied to 
the study of the decay of the inflation in alternative de Sitter invariant vacua |3Cj. 

The Minkowski space-time computation of the decay rate is not suitable for a stage of rapid expansion (as quantified 
above): the rapid expansion of the Universe and the manifest lack of a global time-like Killing vector allow processes 
that would be forbidden by energy conservation in Minkowski space-time. As emphasized in |27ll31| . the lack of energy 
conservation in a rapidly expanding cosmology requires a different approach to study particle decay. The correct decay 
law follows from the relaxation in time of the expectation value of the field out of equilibrium. The relaxation of 
the non-equilibrium expectation value of the field is computed in ref. |27| using the dynamical renormalization group 
(DRG) which allows to extract the decay law directly from the real time equations of motion. The reliability and 
predictive power of the DRG has been tested for a wide range of physical situations including hot and dense plasmas 
in and out of equilibrium j33] • 

The goals of this article: We compute the particle decay of quantum fields minimally coupled to gravity with 
masses M much smaller than the Hubble parameter, which is the relevant case for slow roll inflation. This entails a 
much stronger infrared behavior than for massless particles conformally coupled to gravity. The emergence of infrared 
divergences in quantum processes with gravitons during de Sitter inflation has been the focus of a thorough study 
|33l 1 3 m . As we will see in detail below, similar strong infrared behavior enters in the decay of minimally coupled 
particles with masses M much smaller than the Hubble parameter H. When M « H there is a small parameter 
A ~ M 2 /H 2 which regulates the infrared behavior in de Sitter inflation. We find that a similar parameter A exists 
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in quasi de Sitter slow roll inflation which is a simple function of the slow-roll parameters. A regulates the infrared 
in the self-energy corrections even for massless particles (gravitons). 

We begin by studying the general case of a cubic interaction of scalar particles minimally coupled to gravity allowing 
the decay of one field into two others during de Sitter inflation. The masses of all particles are much smaller than the 
Hubble constant, which leads to a strong infrared behavior in the self-energy loops. We introduce an expansion in 
terms of a small parameter A which regulates the infrared and which in the case of de Sitter inflation is determined 
by the ratio of the mass of the particle in the loop to the the Hubble constant. Long-time divergences associated with 
secular terms in the solutions of the equations of motion are systematically resummed by implementing the DRG 
introduced in refs. [13, E2| and lead to the decay law. We then apply these general results to the case of quasi-de Sitter 
slow roll inflation. We show that in this case a similar small parameter A emerges which is a simple function of slow- 
roll parameters and regulates the infrared behavior even for massless particles. We study the decay of superhorizon 
fluctuations as well as of fluctuations with wavelengths deep inside the horizon. A rather striking aspect is that a 
particle can decay into itself precisely as a consequence of the lack of energy conservation in a rapidly expanding 
cosmology. We then focus on studying the decay of the inflaton quantum fluctuations into their own quanta, namely 
the self-decay of the inflaton fluctuations, discussing the potential implications on the power spectra of primordial 
perturbations and to non-gaussianity. 

Brief summary of results: 

• In the case of de Sitter inflation for particles with mass M <C H a small parameter A ~ M 2 /H 2 regulates 
the infrared. We introduce an expansion in this small parameter A akin to the e expansion in dimension- 
ally regularized critical theories. We obtain the decay laws in a A expansion after implementing the DRG 
resummation. 

• Minimally coupled particles decay faster than those conformally coupled to gravity due to the strong infrared 
behavior both for superhorizon modes as well as for modes with wavelengths well inside the Hubble radius. 

• The decay of short wavelength modes, those inside the horizon during inflation, is enhanced by soft collinear 
bremsstrahlung radiation of superhorizon quanta which becomes the dominant decay channel when the physical 
wave vector obeys, 

k ph (v) < , (1) 

Vv ~ ey 

where rjv, £y are the standard potential slow roll parameters. 

• An expanding cosmology allows processes that are forbidden in Minkowski space-time by energy conservation|27l 
l3lj : in particular, for masses <C H, kinematic thresholds are absent allowing a particle to decay into itself. 
Namely, the self- decay of quantum fluctuations is a feature of an interacting theory in a rapidly expanding 
cosmology. A self-coupling of the inflaton leads to the self-decay of its quantum fluctuations both for modes 
inside as well as outside the Hubble radius. 

• The results obtained in de Sitter directly apply to the self decay of the quantum fluctuations of the inflaton 
during slow roll (quasi de Sitter) expansion. In this case, A is a simple function of the slow roll parameters. For 
superhorizon modes we find that the amplitude of the inflaton quantum fluctuations relaxes as a power law n r 
in conformal time. To lowest order in slow roll, we find T completely determined by slow roll parameters and 
the amplitude of the power spectrum of curvature perturbations A^: 

r= S ^ A \ [l + 0{e v , Vv )\ (2) 



(ey - n v ) 2 

where r\ is conformal time and £y , r\y , ey are the standard slow roll parameters. As a consequence, the growing 
mode which dominates at late time evolves as 

■ (3) 

V 

featuring an anomalous dimension T slowing down the growth of the dominant mode. 

The decay of the inflaton quantum fluctuations with wavelengths deep within the Hubble radius during slow roll 
inflation is enhanced by the infrared behavior associated with the collinear emission of ultrasoft quanta, namely 
bremsstrahlung radiation of superhorizon fluctuations. The decay hastens as the physical wavelength approaches 
the horizon because the phase space for the emission of superhorizon quanta opens up as the wavelength nears 
horizon crossing. 
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• We discuss the implications of these results for scalar and tensor perturbations, and establish a connection with 
previous calculations of non-gaussian correlations. 

The article is organized as follows: In section II we introduce the models, in section III we present the equations of 
motion, describe the approach to obtaining the decay law via the DRG and introduce the A expansion. We consider 
the decay of a scalar field coupled to other scalar fields via a cubic coupling in pure de Sitter inflation. We study 
the decay of super horizon fluctuations as well as of fluctuations with wave vectors deep inside the Hubble radius. 
In section IV we apply the results of section III to the self-decay of the inflaton quantum fluctuations during quasi 
de Sitter slow roll inflation. In section V we discuss the implications of our results for scalar and tensor metric 
perturbations as well as the connection between the quantum decay processes and the emergence of non-gaussian 
correlations. Our conclusions and further discussions are contained in section VI. Two appendices are devoted to the 
calculation of the self-energy kernel in the A expansion for arbitrary wave vector including the order A . 



II. THE MODELS 



We consider a general interacting scalar field theory with cubic couplings in a spatially flat Friedmann-Robertson- 
Walker cosmological space time with scale factor a(t). The cubic couplings are the lowest order non-linearities. Our 
study applies to two different scenarios, i) the inflaton (f> coupled to another scalar field <p, ii) the inflaton field 
self-coupled via a trilinear coupling. We consider the fields to be minimally coupled to gravity. 

In comoving coordinates the action for case i) is given by 

A= [ d 3 x dt a 3 (t)\ ~ (j? - ^J 3 } ~ j- M 2 4> 2 + \ p 2 ~ ^Y^i - \ m 2 tp 2 - g 4> tp 2 + J '(f) <j) + higher nonlinear terms 1 
J 2 2gt 2 2 2a z 2 

(4) 



and for the case ii) 



A = [ d 3 x dt a 3 (t) \ - tj) 2 - - - M 2 <\> 2 + - cj) 3 + J(t) <\> + higher nonlinear terms 1 (5) 

J 2 2gt 2 3 



The linear term in <f> is a counterterm that will be used to cancel the tadpole diagram in the equations of motion. 
The higher nonlinear terms do not affect our results but they are necessary to stabilize the theory. 

In order to study the decay of particles associated with a field we must first obtain the self-energy corrections to the 
equations of motion. We study the decay of inflaton fluctuations up to one-loop order in the coupling either into other 
fields or in self-decay. The calculation of the self-energy correction up to one loop is similar in both cases, the extra 
factor 1/3 in the trilinear coupling in Eq. (JjJJ accounts for the combinatorial factor in the corresponding Feynman 
diagram. Fig. ^ shows the self-energy contributions to the inflaton propagator up to one loop. Fig. displays a 
loop of <p particles and fig. [T\b displays the one-loop self energy for a cubic self-interaction which is the lowest order 
non-linearity around the classical inflaton (expectation value) driving inflation. 

It is clear from these figures that we only need to obtain the self-energy in only one of the cases, since one case is 
obtained from the other by a simple replacement of the masses of the particles that run in the loop. Therefore, we 
will study the self-energy for the case of figure fig. a. 

It is convenient to pass to conformal time rj with drj — dt/a(t) and introduce a conformal rescaling of the fields 

a(t) <p(x, t) = x(x, J]) ; a(t) <p(x, t) = S(x, 77) (6) 
The action Eq. |@J (after discarding surface terms that do not affect the equations of motion) reads: 
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with primes denoting derivatives with respect to conformal time 77, C(rj) = a(t(rj)) is the scale factor as a function of 
77 and 

M 2 x (r 1 )=M 2 C 2 (r 1 )-^ , M 2 {rf) = m 2 C 2 (rj) - ■ (8) 
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FIG. 1: Self energies: Fig. (a) depicts the self energy contribution from a loop of <p particles, Fig. (b) depicts the self-energy 
from the self-interaction of the inflaton. 



In this section we focus on a de Sitter inflationary cosmology, we treat slow-roll and quasi de Sitter inflation in sec. 
IV. For de Sitter space time the scale factor is given by 



o(t) 



r,Ht 



with H the Hubble constant and the conformal time 77 is given by 

V = 



H 



(9) 



(10) 



where rj = — corresponds to the initial time t = 0. 

The Heisenberg equations of motion for the Fourier field modes of wave vector k in the free (<? = 0) theory are given 

by 



xUv) 



k 2 - 



v 1 - 



W+I* 3 - 3(^-4)]*^) = 



where 



9 m 2 
4 ~ IP 



(11) 
(12) 

(13) 



The Heisenberg free field operators can be expanded in terms of the linearly independent solutions of the mode 
equation 



1 



k*--[v 2 -- 

T] 1 



1 



S v (k\rj) = 0. 



(14) 



for v, v respectively. We choose the usual Bunch-Davies initial conditions for the mode functions, namely the usual 
plane waves for wavelengths deep inside the Hubble radius \k r\\ ^> 1. The mode functions S v (q,r]) associated with 
the Bunch-Davies vacuum are given by 



Su{k,ri) = ^i-"-^^H^{krf) . 



(15) 



For wavelengths much smaller than the Hubble radius, ( \krj\ 3> 1) the mode functions with Bunch-Davies vacuum 
initial condition behave as plane waves in Minkowski space time, namely 



Sv (k; V ) '^ie^ 



2A- 



(16) 
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krj 



(17) 



In particular for v = 3/2, ea. (|15f) becomes 

The spatial Fourier transform of the free Heisenberg field operators Xgiv): %( 7 ?) are therefore written as 

Xkiv) = "fc S v {k\rf) + S*(k;r)) 

S s (v) = ^S p (k; V )+^_ s St(k; V ) (18) 

where the Heisenberg operators ag; at and /?■:; j3t obey the usual canonical commutation relations. The Bunch-Davies 
vacuum state |0) is annihilated by a.t,Pz- 

III. EQUATIONS OF MOTION, DYNAMICAL RENORMALIZATION GROUP AND DECAY LAWS: 

THE A EXPANSION IN DE SITTER INFLATION. 

As described in detail in (2^, in a rapidly expanding cosmology the notion of decay 'rate' requires a careful analysis. 
In Minkowski space time the decay rate is obtained from the transition probability per unit time, or alternatively 
from the imaginary part of the space-time Fourier transform of the self-energy evaluated on the particle's mass shell. 
In the transition probability, the square of the energy-momentum delta functions accounts for the space-time volume 
times an overall delta function of energy-momentum conservation: the transition probability divided by this volume is 
the decay rate. In Minkowski space-time, the self-energy is a function of the difference of the space-time coordinates 
due to translational invariance. Hence, a space-time Fourier transform is available, from which the imaginary part 
is extracted. Energy-momentum conservation is of paramount importance to define the decay rate in Minkowski 
space-time, and to determine the kinematic thresholds for particle production and decay. 

In a rapidly expanding cosmology the lack of a global time-like Killing vector prevents energy momentum conser- 
vation, although energy is covariantly conserved. As a result: i) a new definition of the decay 'rate' that does not rely 
on energy-momentum conservation, and a different approach to studying the decay law is necessary ii) since energy 
is no longer conserved, novel processes are possible which are forbidden in Minkowski space-time, therefore we expect 
novel decay channels which are absent in Minkowski space-time. 

In ref.[27j the decay of a particle into massless conformally coupled particles was studied as a test example to present 
the main concepts: the mode functions are the same as in Minkowski space time, this simplified the calculation of the 
self-energy kernel, allowed a systematic study of the reliability of the dynamical renormalization group method and 
a direct comparison to the Minkowski limit. This simple case, however, does not feature several important aspects of 
the more relevant situation of the dynamics of quantum fields which are massless or nearly so but minimally coupled 
to gravity. This latter case features infrared divergences which are not present in the simpler case of conformally 
coupled massless fields [2^, [33[. In this article we study the case in which the inflaton is massive and minimally 
coupled to gravity which is precisely the relevant case for studying the decay of quantum fluctuations during slow roll 
inflation. 

While the main aspects of the dynamical renormalization group method to study decay were introduced in |27| , we 
briefly highlight here the main aspects relevant to this work. 

The method relies on studying the real time relaxation of the expectation value of a field induced by an external 
source in linear response: as the source is switched-off the expectation value relaxes revealing the decay law of the 
amplitude. While an exact solution of the equations of motion is readily available in Minkowski space time because 
the self-energy is a function of the difference of the time coordinates (allowing the use of Fourier-Laplace transforms), 
this is , in general, not the case during inflation. Solving the equation of motion in a perturbative expansion in the 
coupling, secular terms emerge, these terms grow in time when the conformal time r\ — > 0, limiting the validity of the 
perturbative expansion. The dynamical renormalization group precisely allows a systematic resummation of these 
secular terms and the uniform asymptotic expansion provided by the resummation lead to the identification of the 
decay law of the amplitude. Thus, the main steps of the method are the following: 



First, obtain the (retarded) equations of motion for the expectation value of the field in linear response after 
switching-off the source that induces the expectation value. 



• Second, obtain a perturbative expansion of the solution in terms of the coupling. Such perturbative solution 
features secular terms, namely terms that grow in time (when conformal time n — > during inflation) and limit 
the validity of the perturbative expansion. 
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• Third, implement the dynamical renormalization group (DRG) to provide a systematic resummation of these 
secular terms. The solution of the DRG equation gives the decay law of the amplitude of quantum fluctuations. 

We implement these steps in the general case described by the action Eq. Q which couples two fields: <j> and <p 
with masses M and m, respectively and a cubic coupling g<f> ip 2 in exact de Sitter space-time. In section IV we apply 
these results to the case of the self-coupling of inflaton fluctuations. 

We focus on obtaining the decay law for the quantum fluctuations of the field (f> which will be later identified with 
the inflaton field. We derive the equation of motion for the expectation value of the field using the non-equilibrium 
generating functional which involves forward and backward time evolution, typical of a density matrix. Unlike the 
S-matrix case (which is an in-out transition where only forward time evolution is required), the time evolution of an 
expectation value is an initial value problem which requires an in- in matrix element. Real time equations of motion 
obtained from the non-equilibrium generating functional are guaranteed to be retarded. 

It is convenient to write the spatial Fourier transform of the conformally re-scaled field x as follows 

x£fo)=x s (»7) + o£(»j) ; <*£(»?)> =*gfa) ; =o (19) 

where the superscripts ± refer to the forward and backward time branches in the non-equilibrium generating functional 
respectively. The expectation value is the same for both branches since the c-number external source is the same. 
The equation of motion for the expectation value X^{rf) is obtained by requiring (cr^ (j))) = systematically order by 
order in perturbation theory. This is the basis of the tadpole method to obtain the equations of motion, which up to 
0(g 2 ) (one-loop) are given by [see Eq.(37) in ref.[27|] 



2 a 1 f v drf 



r, H 2 



(20) 



where the counterterm J(rj) in the action Q has been used to cancel the tadpole term proportional to (S 2 (x, r/)}, this 
is independent of Xt(rj) and acts as a source term in the equation of motion, v and v are given by Ea. (|13l ). 
The retarded one-loop self-energy kernel /Cp(fc; 7], rj') is given by [27j 



fCp{k;r],ri') 



d 3 q 
(2^ 



Im 



Sp(q,v)St(Q,v')Sp(\q~ k\,f])S*,{\q- k\,rj') 



(21) 



and is depicted in fig. Q}a. The mode functions Sp(q,rj) are given by Eq. (fl5|l . We consider M 2 /H 2 -C 1 and 
m 2 /H 2 <C 1 which for the inflaton case to be studied below corresponds to the slow-roll approximation, and define 



A = -—7 + "777 



3H 2 



H 4 



(22) 



hence A<1. This small parameter A will be related with the slow-roll parameters and plays an important role in 
regulating the infrared behavior in the self-energy. Anticipating a renormalization of the inflaton mass we write 



M 2 = M 2 R + g 2 8 Mf + 0(g 4 



- 2 = 4-s 2 ^ + e% 4 ) 



H 2 



with v\ = 9/4 — Mft/H 2 . The equation of motion up to order g 2 becomes 



XM+9' 



SM 2 



2g 2 
V H 2 



V %lC D {k;r ] ,r 1 ')X % {r 1 ') = Q 
no V 



(23) 



(24) 



In what follows we suppress the subscript R to avoid cluttering the notation, therefore v and the mass must be 
understood as the renormalized ones. A perturbative solution of Ea. l2(J|) is obtained by writing 

X^rf) = X qX (v) + 9 2 X hn (r,) + C( 3 4 ) (25) 

and comparing powers of g leads to a hierarchy of coupled equations: up to second order in g, (one loop order), they 
are 



k -- 



V "4 



v 

4 



(26) 
(27) 
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where the inhomogeneity is given by 

n x {k,ri) = - 



H 2 T] 2 



f v dri 



H 2 Tj 



(28) 



The mass counterterm 5M 2 is fixed by requiring that it cancels the term proportional to X Q k (r])/r) 2 arising from the 
integral in Eq. Ij28(l . The first order correction to the solution is given by 



where Q{k;r),rj') is the retarded Green's function obeying 



dr] 2 



k 2 -- 



v 

4 



Q v {k;r],rf) = S(ri - rj') , Q v (k; 77, 77') = for rj > r\ 



(29) 



(30) 



To compute X 1 ^(77) we first need the kernel ICp(k; 77, 77'). As it will become clear in the following, this kernel involves 
(logarithmic) infrared divergences for v = 3/2, but it is an analytic function of A that features simple poles at A = 0. 
Since in slow roll A<1 we will use the parameter A as a regulator much in the same manner as the e expansion in 
dimensional regularization. We will therefore compute the kernel /Cp(fc; 77, 77') at leading order in A by extracting the 
poles and the logarithmic terms; terms proportional to powers of A give subleading contributions. This is akin to the 
minimal subtraction in dimensional regularization. 



A. Secular terms, DRG and decay law: 

Let g u (k,rj) and fv{k,rj) be two independent solutions of the zeroth order equation lj26|l . the most general solution 

is 

X o,k(v) = H 9u{k, v) + B- k f v (k, 77) (31) 

where and are arbitrary constants. The linear structure of the perturbative series indicates that the perturbative 
solution of the equation of motion has the form 

X % {rj) = A- k g^k^l+g 2 Fi(M) +0(g 4 )} + B % f v (k, rj)[l + g 2 H^r,) + 0(g 4 )} (32) 

The functions Fi(k,rf); Hi(k,rf) are determined by the first order solution Ea. (|29(l and they feature secular terms, 
namely divergent terms in the limit 77 — > 0. Therefore we write, 

F 1 (k;r 1 )=F hs (k;r 1 )+F 1 j{k;r 1 ) ; H x {k;ri) = H hs (k;r]) + #i,/(fc;r/) (33) 

where _Fi iS (fc; 77), H\ tS {k] 77) are secular terms, whereas F\j{k; 77), H\j{k; 77) remain bounded as functions of conformal 
time. The dynamical renormalization group absorbs the secular terms into a renormalization of the amplitudes A^, Bj: 
at a given time scale 77, (wave-function renormalization) |27t l32j | , namely 

A % = A- k (?f) 4(t7) , Z£(v) = l+g 2 z? rk (v)+0(g i ) (34) 
B % = B i (ftZ»<$) , Zf(77) = l+. 9 2 z^(77)+0( 5 4 ) (35) 

The coefficients z^-irf), z^(rj) are chosen to cancel the secular terms in the perturbative solution at 77 — rj order by 

order in the perturbative expansion. Since the scale 77 is arbitrary and the perturbative solution does not depend on 
this scale, the renormalized amplitudes Ar(rf), Bj:(rj) obey the following dynamical renormalization group equation 
to lowest order (Eqs.(55)-(56) of ref . [H H| ) , 



dMV 2 , ^ dF 1>a (k;rj) , ^ 4 

~i)<r 9 ' : " 



9 2 MV) +0(g i ) = (36) 

^- 5 2 ^^#^ + 0(. 9 4 ) = O (37) 
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The solution of these DRG equations is given by 

A^rj) = A u (jh) e?&-'toV- Fl -'te™+°^ , Bffi = B^o) e^A^-H.A^+otf) (3g) 

Setting rj = 77 we obtain the renormalization group improved solution, 

X n ( V ) = A n {rf) g u (k; V )[l + g 2 F.jik^^+Oig^+B^) /^^[l+g 2 H.jik^^+Oig^)] (39) 
A % {r,) = A % {r, ) e » 2 ^..(^)- F M(Mo)] ; B ^ = B _. {m) ^[^(M)-^,^)] . (40) 

The terms in the brackets in Eq. I|39|) are truly perturbatively small at all conformal times. The real part of the 
exponential factors in the complex amplitudes Eq. 14UI) determine the decay law of the amplitude (or growth law in 
the case of instabilities). 



B. The A expansion: 



When the inflaton decays into minimally coupled massless particles, infrared divergences in the self-energy kernel 
are present [29|. These divergences are a hallmark of minimally coupled massless particles namely v — v ~ 3/2 in the 
intermediate state, and are similar to those found for gravitons in de Sitter sr>ace-time|33j. We are instead considering 
the case in which both the inflaton and the decay products are massive with masses M and m respectively and 
(M 2 /H 2 , m 2 /H 2 ) <C 1. The mass of the particles in the loop cutoffs the infrared divergences here, (since m 2 /H 2 ^ 
then P 3/2 and A ^ [Eq. I|22|) ] % 1. As it will become clear in the explicit calculations below, the infrared divergences 
in the self-energy kernel manifest as simple poles at A = 0. Thus, < A <C 1 emerges as an infrared regulator akin to 
the dimensional regularization parameter e = 4 — D in the loop expansion in D-dimensional Minkowski space-time. 

Since A <C 1 for slow roll inflation, we compute the self-energy kernel in an expansion in A keeping the poles at 
A = and the leading logarithms in r/ just like the e expansion in dimensional regularization. (The leading log rj 
terms are the remnant of the infrared divergence regulated by A just like in the e expansion of critical phenomena). 

The details of the calculation of /Cp(fc, 77, 77') in the A-expansion for arbitrary k is presented in the appendix and 
the result for the kernel is given by equations (|B13|I and (|B14|) . We now have all the elements necessary to study the 
decay law. 



C. Superhorizon modes: k — 

We begin by the superhorizon modes and take k = 0. The general solution of the unperturbed mode Ea. (|26[l and 
the retarded Green's function Gv(0, rj, rj 1 ) Eq. (|30H for k = are given by 



X Q M = A (-„)/»+ + B (-r)f- ; & 



± v. 



(41) 
(42) 



We compute the kernel /C p (0, 77, 77') in Appendix A highlighting the most relevant physical processes and displaying 
the origin of infrared divergences as poles at A = and the leading logarithms in r\. /C p (0, 77, ?/) can also be obtained 
in the limit k — > of the k ^ kernel treated in Appendix B [Eq. l|B13Jl ] . 
For k = the self-energy kernel is given by 



alq Im[[S„(q,ri)S*(q,ri')] 2 } , v 



(43) 



The infrared divergences at A = | — v = Q arise from the small momenta behavior of the integrand in Ea. (|43|l . 
Keeping A small but nonzero, we find for the kernel (see appendix A) 



IC 9 (0;r l> r,') = K i (0;ri,r,') + -^. 



1 2 

2A + 3 



i] 2 rj 2 



1-1 

7/j 



V- (44) 
V 77 



where 



1 



77 — 77 



77-77 



8tH (77 - rj'y + (e if) 2 



1 1 

"8^ 2 



1 



1 



T) — if + ier)' rj — rj 1 — it r/ 1 



4 . 

(45) 
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This prescription for the principal part regulates the short distance divergence in the operator product expansion 
with a dimensionless infinitesimal e independent of time. This choice of regularization is consistent with the short- 
distance singularities of the operator product expansion in Minkowski space-time, and leads to a time-independent 
mass renormalization (see also ref. |27|'). 

The two terms in the Eg . (|44l) . namely ICi (0; 77, 77') and the term in braces have very different origin. K-i (0; 77, 77') 
accounts for the large loop momentum contribution qrj, qrj' 3> 1 where the behavior of the mode functions is the same 
as for conformally coupled massless fields, in particular, the short distance (ultraviolet) divergence is present only in 
this term. The terms in the braces account for the strong infrared behavior of superhorizon wavelengths reflected by 
the pole at A = and the logarithms in 77. This calculation exhibits clearly the origin of the different contributions. 



-(f ~ 



/c = 



<f ~ 

FIG. 2: Decay of the field <j> (solid line) for k — into superhorizon modes of the field <p 

(dashed lines) 



It remains to integrate over 77' in the second term in Eq. (|28l) with the kernel given by Eg . (|44|) . The integral 
involving K. 1 (0; 77, 77') in Eg . (|44l) was given in ref. 27]. The integrals over 77' for the second term (between braces) in 
Eo . (|44|l can be done easily by expanding the ln[l — 77/77'] in a power series in 77/77' and integrating term by term. The 
result of the integral over 77' in Eq. I|28[l is of the form 



7J- /" ^ MO; V, rf) Xojtf) = A (~ v f+ 2+ + B {-rjf- ^ + F[ V , m ] 



(46) 



where ^[77, 770] refers to the contribution of the lower integration limit and does not produce secular terms in X^(r\). 
Integrating over 77' in Eo. (|46|l yields, 



1 



(2nH) 



lne + 7 + 7/>(- - v) 



3{itH) 



,2 _ 9 V 2 A 

4 v 



(^ + |) 2 



2 " 2 ^ " 

_ (47) 

where we used Eg . (1441) . 7 is the Euler-Mascheroni constant, and the contribution from K. 1 (0, 77, 77') [23. a_ follows 
from a + by changing v —* —v while v is unchanged. Introducing the symmetric and antisymmetric combinations, 



a s = - [a + + ol- 



a a = ^ i a + _ a - 



the symmetric term, 



[A (-77)^+ + B {-vf-} 



(48) 



(49) 



is identified with a contribution to mass renormalization and is cancelled by the counterterm SM^ including the 
logarithmic ultraviolet divergence lne [e the short distance regulator Eq. (|45|l ] . Eos. (|47|) and l|48|l yield, 



1 



8vriJ 2 



tan[7r v\ 



(50) 



The unit term in the bracket arises from the contribution of K, 1 (0, 77, 77'). After cancelling the term given by Eq. I|49|) 
with a proper choice of the mass counterterm, and taking into account that F[rj, rj ] does not contribute to the secular 
terms, we find the solution of the equation of motion up to 0(g 2 ) 



X g (r?) = X 5(77) 



1 + r In h non — secular terms 

?7o 



(51) 
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with 



16tt B 2 v 



tan[7r u] 



g 



16tt B 2 v 



tan[7r u] 



AH 2 



M 2 



M 2 
Jl 2 



(52) 



The dynamical renormalization group resummation exponentiates the secular terms in Ea. (|51|l [see Eas. (|39|) - (I40|) ] and 
leads to the improved solution, 




A{m) (-V) P+ [1 + 0{g 2 )] + B(r, ) {- V f- [1 + 0(g 2 



(53) 



The first term inside the square bracket in Eq. I|52[l (namely the unit term) corresponds to the case in which the 
inflaton decays into massless particles conformally coupled to gravity [2 ^ |.[2 ^ |. 

The calculation leading to ea. (|52|l is valid for v — > | (namely, m <C H) and we keep v as well as M arbitrary. We 
can analytically continue the formula l|52|l to H < M and then take the m <C H <C M limit. In this limit T becomes 
the decay rate of a particle with mass M into massless particles in Minkowski space-time: 

lira H T = T Mmk = g 2 /{16nM) , 



as it must be. 



D. Modes inside the horizon during inflation: \kn\ S> 1 



The kernel K,p(k; 77, 77') for arbitrary k has been computed in the appendix in leading order in the A expansion 
and up to leading logarithms. It is given by equations 1|B13(1 and (|B14I) . Obtaining the perturbative solution and 
extracting the secular terms leading to the decay law for arbitrary k is an extremely difficult task which requires a 
full numerical study. However, explicit expressions can be derived for wavelengths deep inside the Hubble radius all 
throughout inflation, namely \krf\, \krf\ >> 1. In the appendix we show that in the short wavelength limit the kernel 
simplifies to the following expression [Eq. i|B16(l ] 



ICp(k, 77,77') = K.i (fc, 17,77') - 



Air 2 krjrf 



-\ sinfc(7?- 77') 



+ C - In fc(77 - ?/) + In fc 2 77?;' - Ci[2fc(?7 - 77')] 



+ 003/0(77 — 77') — + Si[2fc(?7 — 77')] 



(54) 



where C = In 2 + 7 — 2. Again, the first term in this expression is the self-energy kernel for conformally coupled 
massless fields in the loop, 



ICi(k,T],r]') 



7r cos fc(?7 — 77') V 

8tt 2 



77-77' 



(55) 



The principal part is defined by Eq. I|45|l . A close examination of the steps leading to this expression in the appendix 
shows that this contribution originates solely from the high loop momenta qrj, qrf 3> 1 for which the mode functions 
coincide with those of massless conformally coupled fields. The second term within brackets which features the 1/A 
and the logarithms originate in the process of emission of supcrhorizon modes. Two regions in the integral over the 
loop momentum q give rise to these contributions: q I/77 and \k — q\ <^ 1/7/, corresponding to the case when 
either line in the loop transfers very small momentum (superhorizon modes). These processes can be described as 
bremsstrahlung radiation of superhorizon quanta and are depicted in fig. |3J 

This process is analogous to the generation of Hawking radiation from black holes. In the case of Hawking radiation, 
a pair is created from the vacuum, a particle falling inside the horizon and the other one being emitted outside. In 
the present case a particle inside the Hubble radius decays into a pair: a particle goes outside the Hubble radius 
and the other inside. Analysis of the phase space integration carried out in the appendix reveals that the emitted 
superhorizon quanta are almost collinear with the (large) external momentum k. 

In summary: the processes which yield the leading infrared behavior responsible for the term 1/A and the leading 
logarithm in the kernel Eq. I)54|l correspond to collinear bremsstrahlung radiation of superhorizon quanta. 
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-» 7t k — q 

A: / k ' 
< + < 

\ s 
s \ 

FIG. 3: Infrared contributions to <f> — > The external particle has a wavelength deep inside the horizon but one of the 

intermediate lines has superhorizon wavelengths. This process is identified as bremsstrahlung radiation of superhorizon quanta. 



In the limit \kr)\ >> 1 the zeroth order equation of motion is 

^(l) + ^o|(l) 



(56) 



whose solutions are simple plane waves describing short wavelength modes deep inside the Hubble radius, 

XoM=A j ;e-^+B % e^ (57) 



The equation of motion for the first order perturbation is given by 

X'l j .(r 1 ) + k 2 X lj :( v )=TZ 1 (k,f 1 ), 



(58) 



where the inhomogeneity is given by Eq. (|28J) with X Q ^(rj) given by Ea. (|57|l . The calculation of lZi(k, 77) for general 
k is very complicated, but the leading order terms in the limit \krj\ ^> 1 can be extracted systematically. There are 
two distinct contributions : i) the first term in Eq. (|54|l , namely the short wavelength modes which yield the kernel of 
conformally massless fields /Ci (fc, 77, rj'), and ii) the superhorizon modes yielding the second term in Ea. H54|l with the 
1/A and the leading logarithms. 

The first term contains a short distance divergence proportional to X Q ^(77) [lnej/r? 2 where e is defined in Eq. I|45|l . 
This term is canceled by the proper choice of the mass counterterm in the inhomogeneity Eq. I|28fl . The second 
term does not yield a mass renormalization to leading order in A. After a proper choice of the mass renormalization 
counterterm, we find 



7T ?7 

i In — 

2 Vo. 



4w 2 kH 2 r] 3 I k 



1 i / 1 \ . 7T 

— + ln(-fcr/) - 1 - 



- + 1 In — 

2 m 



B e ikri 
k 



1 1 / , s ■ 7r 

— + ln(-fc?7) + 1 - 



(59) 



where we have displayed separately the contributions from the high momentum modes yielding the first term and 
arising from K, 1 (k, rj, 77'), and those from the superhorizon modes yielding the second term which feature the hallmark 
1/A and logarithms. We have kept the leading order terms in the real and imaginary parts inside the brackets 
neglecting terms suppressed by higher powers of l/|fcr/| << 1. A noteworthy feature of the contribution of the 
superhorizon modes is the extra factor 1/krj. The origin of the extra factor 1/k can be traced from the phase space 
angular integration Eq. I|B2(I . The integrals yielding JC 1 (fc; 77, 77') are dominated by momenta q > k which compensate 
the factor 1/k. The integration over the superhorizon modes cannot compensate the 1/k for large k corresponding 
to modes well within the horizon. Hence, the extra factor 1/k is a consequence of the small phase space available for 
the coupling between high and small momentum modes. For dimensional reasons this extra factor k appears with 
an extra factor 77 which is the only other scale in the integrals. In summary: the contribution of the superhorizon 
modes yields a strong infrared behavior which is regulated by A and is suppressed by phase space by an extra power 
of 1/|A;77| = HC{ri)/k <C 1 in the limit |fc77| ^ 1, (i. e., for wavelengths much smaller than the Hubble radius during 
inflation). 
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The solution of the equation of motion (|58|l is found by using the retarded Green's function in the short wavelength 
limit, which is given by 



Gu(k;r],r]') = ^ sin k (77 - rf) 6(77 - 7/) 



Therefore, the first order correction becomes 



sin k(rj — ?/) IZi {k.-q') drj' . 



(60) 



(61) 



The final integral with the retarded Green's function as in Eq. (|61|) can now be performed extracting again the leading 
order terms for \krj\ 3> 1. We find 



1 



1 Z2irkH 
secular 



[C(v) - CM] 



1 



167T 2 fc 2 



f -'(,/) (i + ln( •/.•//) I - r'(, 7ll ) (--: 1„( -/.•//,,) 



(62) 



where C{rf) = —1/Hr) is the scale factor and we have omitted purely imaginary secular terms since the dynamical 
renormalization group exponentiates them [Eci. (|38|) ] to a (time dependent) phase. The terms displayed in Ea. l|62|) are 
truly secular, since log \rj\ grows by about 60 during inflation. From the dynamical renormalization group resummation 
Eq. (|38[) we find the following improved solution of the equations of motion 

X%in) = A % {ri) e-^ [1 + <D(g 2 )] + B^rj) [l + 0(g 2 )] 

B % {r}) = B s { Va ) e- [rik '^- r< - k ^>)] e ~ 1 ^) (63) 
where the real phase £ (fc, 77) is not relevant for the decay rate, and the decay law of the amplitudes is given by 



T(\kj]\ > 1) = 



H 



327riJ 2 hphtn) 



1 



2^ fl ln k ph ( V ) 
irk ph (rj) \A H 



k P h{r]) 



(64) 



Modes with k deep within the horizon satisfies k p h(r])/H ^> 1. The unit term in the bracket in Eq. 164|) corresponds 
to the particles in the loop being massless and conformally coupled to gravity. The second term which features 
the factors 1/ A and the logarithmic term arise from the emission of superhorizon quanta. We see that the infrared 
regularization provided by A yields a finite result for the decay law. Furthermore, the factors associated with the 
infrared processes are suppressed by an extra power of H/k p h(j]) <C 1. This suppression is a consequence of the small 
phase space available for the coupling between high and small momentum modes as can be seen directly from Eq. 
(lB2l . 

The contributions to the decay law Ea.l|64H from the emission of superhorizon quanta become larger the closer is 
the wavelength to horizon crossing. For H/k p h(f]) 1 they can even dominate T for sufficiently small A. 

An important aspect of the decay law, either for modes inside or outside the Hubble radius, is that there are 
no kinematic thresholds. This is a consequence of the inflationary expansion [(M, m) H] and the lack of energy 
conservation. In Minkowski space-time, energy-momentum conservation leads to kinematic thresholds, in particular 
a massive particle cannot decay in its own quanta. However, in an inflationary cosmology this process is allowed, 
namely a particle can decay into itself. 



IV. SELF-DECAY OF INFLATON QUANTUM FLUCTUATIONS DURING SLOW-ROLL INFLATION. 

Fluctuations of the inflaton in exact de Sitter inflation do not seed density perturbations: scalar metric perturbations 
couple to the inflaton fluctuations through the time derivative of the inflaton expectation value. Therefore, the relevant 
case for density perturbations is quasi de Sitter inflation, in particular slow-roll inflation 0, |2^| , which serves as the 
basis of CMB data analysis 16]. Our ultimate goal is to understand how quantum effects from interactions, such as 
the decay of fluctuations, can affect the power spectrum of scalar and tensor metric perturbations. During slow roll, 
the scalar field fluctuations are sources for scalar metric fluctuations, therefore quantum effects as studied here can 
produce novel signatures on the power spectrum. While a complete gauge invariant description is ultimately required 
to treat this issue, we focus here on the decay of the inflaton quantum fluctuations, which in longitudinal gauge are 
directly related to the scalar metric fluctuations 0, |3^| . 
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Therefore, we apply our results above to the quantum self decay of the inflaton fluctuations. We consider only one 
scalar field, namely the inflaton whose Lagrangian density is given by 

C=- d^d^ - V[$] (65) 

We write 

= #o(t) + (66) 

where t is the cosmic time, §o(t) is the expectation value of the inflaton field which drives the FRW background metric 
and <f>(x,t) describes the inflaton quantum fluctuations. Expanding around the expectation value, the Lagrangian 
density for the quantum fluctuations reads 

8C[<j>; *„] = r (d"<P - V"[<S> ] <p 2 ) + f <t? + higher order terms , (67) 

where the primes applied to the potential V[$] stand for derivatives with respect to the argument (not to be confused 
with derivatives with respect to conformal time) and we have used the equation of motion for $o(0 which in cosmic 
time is given by 

0o + 3iT$o(*) + V r/ [*o] = O. (68) 
We have kept the lowest order term in the non-linearity and defined the (dimensionful) coupling constant 

g=lv'"[^ ] . (69) 

As it is clear from the study in the previous section, the perturbative treatment of the non-linearity will be reliable 
provided g/H <C 1. 

In the slow roll approximation the equation of motion simplifies to 

3 if $o(f) + V'[*o] =0- (70) 
The slow-roll parameters relevant to our discussion are the following (either in terms of H (Hubble) or V (potential)) 0, 

m 

F -2M 2 ( H '\ 2 f _M pl (V'[^ ]V_ 

1ff = 2M^ , r, v = M 2 Pl ^^= m + e v (72) 

£, H - 4 M Pl H2 , £, v = M Pl v 2 [$> ] ' 

Here, Mj,, = l/[8nG] = m 2 pl /(8n) and M Pl = 2.4 10 18 GeV. Slow roll implies (e v , r) V , £y, e H , rj H , £ H ) < 1, and 
£y j £p are formally of second order in slow-roll. 

In terms of the slow roll parameters the Friedmann equation reads 



E 2 = 



V 



3M pi 



l + JL + O^evw) (74) 



and the effective mass of the inflaton quantum fluctuations is given by 

M 2 = V"[<S>o}=3 H 2 , lv + 0(e v rj v ) . (75) 
During slow roll inflation the scale factor is quasi de Sitter and to lowest order in slow roll: 

^) = -^^ = -^(l + £ .) + 0(4) (76) 

We are now in a position to apply the results obtained in the previous sections to study the self decay of the quantum 
fluctuations of the inflaton via the cubic coupling. The decay process <j> — > <j> <fi is depicted in fig. The self-energy 
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has the same form as for an interaction <j> ip 2 analyzed in the previous section, the only difference is that in the self 
energy we have v instead of v. 

Within slow roll, the linearized equation of motion for the quantum fluctuations of the inflaton <j> is precisely given 
by 



xUv) 



k 2 + M 2 C 2 {r/) 



C{v) 



xM = o , 



(77) 



with the quasi de Sitter scale factor C(rf) given by Eg. ((76(1 and M 2 given by Ea. (|75|l . Let us compute the expression 
within brackets in Eq. 1(77(1 to first order in slow roll. We have in cosmic time: 



= a 2 (t) 



H + 2W 



And Eqs. (JZOJ, dJ and O yield, 

Eqs.lO, CHI and (2J thus imply, 

C"(V) 



and 



M C 2 {r/) 



c"(v) 



JL = -e v + 0(e 2 v ) 



= — [2 + 3 ey + higher orders in slow roll] 



[2 + 3 (ey — T)v) + higher orders in slow roll] 



(78) 



(79) 



(80) 



(81) 



where we also used Eq. 1(75(1 . In summary, during slow roll the quantum fluctuations of the inflaton behave in a 
similar way as in pure de Sitter space-time [sec. II and II] but with the index v of the mode functions given by 



ey — rjy + higher orders in slow roll 



(82) 



where 



C(rj) - if 

In this case A expresses in terms of the slow roll parameters as 



higher orders in slow roll 



where 



3 1 

A = — — v ^ rjy — ey + higher orders in slow roll = -(ng — 1) — tit 



ns — 1 — 6 ey + 2 r/y , ut = — 2 ey = — r/ 



(83) 



(84) 



are the scalar and tensor spectral indices, respectively and r is the tensor to scalar ratio. Thus, we can apply the 
results obtained in the previous sections to this case in which the inflaton quantum fluctuations decay into themselves 
via the trilinear coupling. To leading order in slow roll this is done simply by setting v — v — 3/2 — (rjy — ey) in the 
results previously obtained. The slow roll parameters remain constant to leading order in slow roll. 

Superhorizon modes: k = 

For superhorizon modes the results Eas. 1)52(1 and 1(53(1 show that the amplitude of the quantum fluctuations decay 
as i] — > with the power law given by Ea. 1(53(1 with /3 + = 2 + ey — r/y, (3- = — 1 + r/y — ey and 

n 2 



n H (r/y — ey) 



[l + 0(ey,r/y)} 



(85) 
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FIG. 4: Self-decay of quantum fluctuations of the inflaton. All lines correspond to the field (j>, i-e, the quantum fluctuations of 
the inflaton. 



The effective dimensionless coupling g/H [see Ea. ()69|l ] is related to the scale of inflation and the slow roll parameters: 

g_ = V"'[<i> ] Jy_ _3_H_ 

H 2 H 2 Mpi 1 ' 

To lowest order in slow roll the power spectrum of curvature perturbations A 2 ^ is given by Q 

H 2 1 V 3 

A ^ = 8 ir 2 M 2 pl e v = 12 tt 2 M e pl ' (8?) 

This allows to relate the effective dimensionless coupling g/H to quantities that are observable from CMB data: 

j I =3n{; v (A 2 z yil + 0( W ,e v )} . (88) 

Then, we can write T [Ea. (|85|l ] completely in terms of slow roll parameters and the power spectrum of curvature 
perturbations, 

r= 8#A» x 

and in terms of the scalar index and the tensor/scalar rate we have using Eq. 184(1 . 

r = 32 {n fl f+ 1)2 t 1 + 0(e v ,riv)} ■ (90) 

Of particular importance is the behavior of the growing mode of the inflaton which gives the dominant term for rj — > 
after the wavelengths of the fluctuations cross the horizon. ^From Eqs. I(41jl . (|53|l and 182(1 . it is given by 

/-+ r = 5 . (91) 

We see that the decay constant F acts as an anomalous scaling dimension for the growing mode of the supcrhorizon 
fluctuations of the inflaton. The decay rate T slows down the growth of the dominant mode for 77 — > and fastens 
the decrease of the subdominant modes. 

Modes inside the Hubble radius during slow roll inflation |fcry| ^> 1: 

The decay law of the quantum fluctuations of the inflaton with wavelengths deep inside the Hubble radius \kr)\ = 
k p h{rj)/H 3> 1 is given by Eas.l(63|) with 



.,2 



r( kn » 1) = "' U 



32irH 2 k ph {r]) 



2H f I t ln k ph ( V ) 



t k ph (rj) \rj V -e v H 



l + 6y 



C(rj) wy Hrj 

(92) 



17 



The term that is inversely proportional to the slow-roll parameters and the logarithm of k p h (rj) /H are a consequence 
of the almost collinear emission of superhorizon quanta. This process corresponds to one of the decay lines in fig. 
0] carrying small momentum with superhorizon wavelengths as explained in the previous section. This process is 
identified as the emission of bremsstrahlung radiation of superhorizon quanta. 

Which of the terms in the bracket in Eq. (|92|l dominates depends not only on the numerical value of the slow roll 
parameters but also on how close the physical wave vector is to horizon crossing. For very short wavelength modes, 
namely 

2 H 

k ph (v) > (93) 

rjv - ev 

the second term is negligible and the result is similar to the inflaton decaying into massless conformally coupled 
particles. This is of course due to the fact that in this kinematic region the modes in the internal propagators are 
simply plane waves with Bunch-Davies initial conditions. On the contrary, as k p h(rj) approaches horizon crossing, 

k ph (v) < — — , (94) 
rjv - ev 

the emission of ultrasoft collinear quanta, namely superhorizon bremsstrahlung radiation, becomes the dominant 
decay channel and the second term dominates. This crossover phenomenom can be interpreted as the phase space 
for collinear emission opening up near horizon crossing. This is because the phase space factor 1/fc in eq. (|B2ll is 
effectively l/|fc?7| = H/kphir/) by dimensional reasons. 

Thus, in slow roll, there is a wide region of physical momenta 

J ff«fc P / l W« — ^— . (95) 
Vv ~ ey 

for which the first term in the bracket in Ea. 1)92(1 can be neglected and the leading slow-roll result Ea. 188(1 for g 2 /H 2 
can be used: 

r(NI»l) = ^ (^t^) (1 + 2 e v ){—^—+\og[k 77(1 - e v )]} . (96) 
16 V k V J I Vv - ev J 

V. IMPLICATIONS FOR THE DECAY OF SCALAR AND TENSOR PERTURBATIONS AND 

NON-GAUSSIANITY 

While we have focused on the decay of the inflaton quantum fluctuations during slow roll we can extrapolate our 
results to see how our findings may provide corrections to the power spectra of scalar and tensor perturbations. 

A. Curvature perturbations: 

For scalar perturbations, the action for the gauge invariant perturbation 

u k = C{rj) <j) k + % Vfc , (97) 
rl 



has a simple form at quadratic order |l2j | and obeys the equation of motion 



c 2 - — 

z 



u k = ; z = . (98) 
rl 



Here <£> is the inflaton expectation value, 4>k the inflaton fluctuation and ipk are the spatial curvature perturbations. 
The gauge invariant perturbation u k is related to the curvature perturbation on comoving hypersurfaces 1Z as [6l ll2Ll2^| 
Uk = — z IZk- Therefore, the power spectrum of the curvature perturbation A 2 ^ is directly related to the corresponding 
spectrum of the gauge invariant perturbation Uk |rlll2ll2^| 
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For superhorizon modes, k 2 <C z"/z, the only relevant contribution is the growing mode Uk — Ak z. Therefore, well 
after horizon crossing the power spectrum of the curvature perturbation is constant in time and given by 



During slow roll H HI 



— = 2 C 2 ( V ) H 2 



1 + e v - - Vh 



1 + ^{W + ev) + 0(T]y,ey : r] V e v ,t v ) 



where we used Eas. lfT^ and l(7B|l and which implies 

3 

v = - + ev + Vv + 0(rjy, e v , r] V e v , £ v ) 
Therefore, the small parameter A for the gauge invariant perturbation is given by 

A = -ev - Vv + 0(r) v ,e v ,i]vev,€v) 

Without interactions among fluctuations (no decay of fluctuations), the growing mode for superhorizon wavelengths 
to lowest order in slow roll is given by 



(100) 



(101) 



(102) 



(103) 



fa) = A k rT (1+£v+7?v) 



(104) 



When a cubic interaction for the perturbation u^irj) is introduced, the results of the previous sections imply that an 
anomalous scaling dimension, namely the decay rate T will appear in Eq. (|104fl . i. e. 

u k {7 1 )=A k r 1 - (1+ " v+m/ - T) , 

as in Eq. Ij91|l for the inflaton fluctuations. Unless z acquires the same anomalous dimensions, the amplitude of the 
power spectrum [Eq. (|99|l ] will depend on time. 

In order to study the decay of curvature perturbations the next step in this program is to obtain the cubic vertex 
for the variable u and to compute the one-loop self energy. This will ensure the gauge invariance of the results. The 
gauge invariant formulation of ref.|l2| has to be extended to higher order of perturbations, as for example in ref.|37|. 
or alternatively, we can work in a fixed gauge. The cubic interaction vertex for three scalars has been computed in 
ref.|2l|. In particular, the contribution to the self-energy from ghost loops must be included. 

The computation of the self-energy corrections will follow the same lines presented in the previous sections with the 
extra feature of momentum dependent vertices and ghost-loops. The infrared behavior of the loops will be regulated 
by A given by Eci. ((TU5|) . 

In order to obtain corrections for the power spectrum of curvature perturbations, the secular terms arising from the 
equations of motion will have to be resummed by DRG for the whole range of physical momenta, until and beyond 
horizon crossing, and establish the behavior of the growing mode. At least two sources of corrections to the index of 
the power spectra may be expected: 

i) from the amplitude of the growing mode Ak , which is obtained by matching the solutions deep within the Hubble 
radius and well after horizon crossing, 

ii) from the anomalous scaling dimension of the growing mode, namely the decay rate F. 

Furthermore, in order to extract the power spectrum, the corrections to the dynamics of the zero mode of the 
inflaton <j> from the coupling of the zero mode to the fluctuations must be obtained. This study will indicate whether 
the variable z also acquires an anomalous dimension and if so, whether the ratio u k /z for the growing mode is time 
independent when the decay process is accounted for. This study is in progress. 



B. Gravitational waves: 

For tensor modes, i.e, gravitational waves the action at the quadratic level is simple |T^ and the field operators for 
the (gauge invariant) gravitational waves are expanded in terms of the mode functions Vk which satisfy 



rill 

h c 



v k = o, 



(105) 
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where %r- in given by Eq. (gj. We then have 

C" _v 2 ~\ 



C T) 



2 



higher orders in slow roll , (106) 



3 

v = — + ey + higher orders in slow roll , A = — ey + higher orders in slow roll . 

Thus, during quasi de Sitter slow roll, the mode functions for gravitational waves are Hankel functions with an index 
v slightly different from 3/2. 



The three graviton vertex was computed in refs. 21, 33] . The Born scattering amplitude for three gravitons in a 
de Sitter space-time features both infrared as well as secular divergences^. These infrared divergences precisely 
arise because the mode functions for gravitons in de Sitter space-time are Hankel functions with index v = 3/2. The 
remaining long time (ry — > 0) divergences can be consider as secular terms akin to those found above. 

During quasi de Sitter slow roll inflation, the parameter A can regulate the infrared divergences found in refs. 
|33l |34|, and the DRG implemented here will provide a resummation of the secular long time divergences found in 
refs. |33j and Q. 

While the three graviton scattering amplitude has been computed in the Born approximation, the full self-energy 
for gravitons has not yet been obtained. In order to understand the decay law for gravitons, the program presented 
in this article must be implemented. Such program for gravitons, as in the case of the scalar perturbations discussed 
above, may require including the contribution from ghost loops to the graviton self-energy. Clearly, carrying this 
program in either case is a task beyond the scope of this article. 

C. Estimates of the cubic coupling and the decay rate from the WMAP data 

In the slow roll approximation the slow-roll parameters themselves are slowly varying functions of time, in particular 



Vv „ /k — 2Mp; g 3/2, 
— = 2e v r, v - V2ey — + 0(g ) . (107) 



The slow roll approximation entails that fjy / H <C 1 which in turn requires 

£ « • (108) 

WMAP gives the following value for A^ 

A 2 n ~ 2.2 x 10~ 9 (109) 

which when combined with Eq.|(S7J yields the following estimate on the scale of inflation 

H - 10 15 y/e^ [Gev] (110) 

The WMAP data on dn s /d\nk suggests that £ v ~ 0.028 ± 0.015 which combined with Eq. JHSJ leads to the 
following estimate on the cubic coupling, 

-|~1.3xl0- 5 (111) 
H 

which places the strength of the dimensionless coupling within the validity of the slow roll approximation and pertur- 
bation theory. These results in turn lead to the following estimate for the rate T [Eq. I|89|) ]. 

3 x 10" 8 > T > 3.6 x 10" 9 . (112) 

This gives in cosmic time for a typical value H ~ 10 14 GeV 

10 7 GeV >T ds = HT> 10 6 GeV . (113) 
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D. Connection with non-gaussianity: 



Non-gaussianity of the spectrum of fluctuations is associated with three (and higher) point correlation functions. 
An early assessment of non-gaussian features of temperature fluctuations in an interacting field theory was given in 
ref.^^]. In ref.0| the simplest inflationary potential with a cubic self- interaction for the inflaton field was proposed as 
a prototype theory to study possible departures from gaussianity. The three point correlation function of a scalar field 
in a theory with cubic interaction as well as the four point correlation function in a theory with quartic interaction 
were calculated in ref . jU ■ 

The long time (rj — > 0) behavior of the equal time three point correlation function in the scalar field theory defined 
by Eq.JSJ for M = (and hence v = §), is given by [3 



(X(k, r,) X % V) Xi-k - q, r,)) = ^ C 3 (n) 9 2 ' " ' ' /: ! 

k q \k + q\ 



(114) 



where 



F(k,q\rj) = k 3 + q 3 + \k + q\ S {\n{k T il)+l]-{k 2 + q 2 + \k + qf)k T + kq\k + q\ ; k T = k + q+\k + q\ (115) 

A diagrammatic interpretation of the equal time expectation value Eq. (|114|) is depicted in fig. [SJ which illustrates the 
similarity with the decay process depicted in fig. 



■k 



q,r] 



k,r] 




q,r] 



FIG. 5: Equal time three point function (x(k,ri) x(<f) 7 ?) x( 
the vertex is integrated. 



q, rj)) in the Born approximation. The time coordinate rj\ of 



Furthermore, the logarithmic secular term in Ea. l|115|l indicates that the three point function features secular 
divergences even at the tree level. It is argued in rcfs. dUil that - \n[k T r}\ - 60 which is the number of 
c-folds from the time when fluctuations of wavenumber k^ first crossed the horizon till the end of inflation. However, 
such infrared logarithms arc secular terms and have precisely the same origin as in the self-energy kernel and in 
the inflaton fluctuations discussed here (sees. IIIC and HID). The same holds for the infrared logarithms in the three 
graviton scattering vertex j^J Ell • 

In particular, the self-energy computation corresponds to a further integration over the loop momentum q. In a 
fairly loose manner, the self-energy is basically the square of the three point correlation function integrated over the 
loop momentum. This is akin to the unitarity relation between the imaginary part of the forward scattering amplitude 
and the square of the transition amplitude in S-matrix theory. 

In summary, the interaction between the fluctuations gives rise to non-gaussian correlations which are determined 
by the three point function which is precisely related to the self-energy and the decay of the quantum fluctuations. 
Therefore, the decay of the quantum fluctuations of the scalar field will also lead to non-gaussian correlations and 
non-gaussianity in the power spectrum. 

The direct relationship between the self-energy, decay and non-gaussian features of the power spectrum will be the 
subject of further study. 
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VI. CONCLUSIONS AND FURTHER QUESTIONS 



In this article we have studied particle decay of fields minimally coupled to gravity in the case when the mass of the 
fields is -C H during inflation. Unlike the decay into massless fields conformally coupled to gravity, this case features 
a strong infrared behavior which leads to novel results. 

We have implemented the dynamical renormalization group resummation program introduced in ref. |27j combined 
with an expansion in a small parameter A which regulates the infrared. 

In the case of exact de Sitter inflation, A is a constant equal to the ratio of the mass of the decay products to the 
Hubble constant, while in slow roll inflation A is a simple function of slow roll parameters. The expansion in A is 
akin to the s expansion in critical phenomena in dimensional regularization. The dynamical renormalization group 
provides a resummation of the long-time secular divergences which determine the decay law of quantum fluctuations. 

The lack of energy conservation in an expanding cosmology leads to the lack of kinematic thresholds for particle 
decay. In particular, this possibility leads to the self-decay of quantum fluctuations whenever a self-interaction is 
present. 

We have studied the decay of a particle for a cubic selfcoupled scalar field in de Sitter space-time and applied the 
results to the self-decay of the inflaton quantum fluctuations during quasi de Sitter, slow roll inflation. We focused 
on extracting the decay law both for wavelengths well inside and well outside the Hubble radius. In both cases the 
strong infrared behavior enhances the decay. 

The decay of fluctuations with wavelengths much smaller than the Hubble radius is enhanced by the collinear 
emission of ultrasoft quanta, this process is identified as bremsstrahlung radiation of superhorizon quanta. As the 
physical wavelength approaches the horizon, the phase space for this process opens up becoming the dominant decay 
channel for short wavelength modes in the region 

H < k ph (r)) < H . (116) 

The decay of short wavelength modes hastens as the physical wavelength approaches the horizon as a consequence of 
the opening up of the phase space. 

Superhorizon quantum fluctuations decay as a power law ~ rf in conformal time, where T is determined by the 
following combination of the slow roll parameters and the amplitude of curvature perturbations 



32 f 2 A 2 

r = („ rrxna t 1 + °kv>nv)\ (117) 

\n s i -f- .) 



This decay law entails that the growing mode for superhorizon wavelengths evolves as rj Vv ~ ev+T /rj hence T provides 
an anomalous scaling dimension slowing down the growing mode for — >• 0. 

The recent WMAP data indicate that 3. x 10~ 8 > T > 3.6 x 10~ 9 . This corresponds to a decay rate in cosmic 
time 10 7 GeV > I^s = H T > 10 6 GeV. Although these values may seem small, it must be noticed that the decay is a 
secular, namely cumulative effect. 

We discussed some potential applications and implications for primordial scalar and tensor perturbations as well 
as the relationship between the decay processes studied in this article and the generation of non-gaussian features in 
the primordial power spectrum. 

The results of our study bring about several questions: 

• The generation of superhorizon fluctuations during inflation is usually referred to as 'acausal'. However, we 
have found that fluctuation modes deep inside the horizon decay into superhorizon modes, therefore there is 
a coupling between modes inside and outside the horizon. The phase space for this process opens up as the 
physical wavelength approaches the horizon. It is natural to conjecture that this process that couples modes 
inside and outside the horizon with a coupling that effectively depends on the wave vector, will ultimately lead 
to distortions in the power spectrum. This distortion will necessarily be small in slow-roll since the coupling is 
of the order of the slow roll parameters, but it may compete with the running of the spectral index from the 
non-interacting theory which is itself of quadratic order in slow roll. 

• In the non-interacting theory, the equation of motion for the gauge invariant Newtonian potential (equal to 
curvature perturbation) features a constant of motion for superhorizon wavelengths [TH IT^ . This is used to 
estimate the spectrum of density perturbations in inflationary universe models. It is conceivable that this 
conservation law will no longer hold in higher order in slow roll when interactions are included. We expect this 
to be the case for two reasons: the coupling between modes inside and outside the horizon as well as the decay 
of superhorizon modes. Clearly the violation of the conservation law, if present, will be small in slow-roll, but 
this non-conservation may also lead to distortions in the power spectrum. 
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• While we have focused on the decay process during inflation, our results, in particular the decay of superhorizon 
fluctuations and the coupling between modes inside and outside the Hubble radius, raise the possibility of similar 
processes being available during the radiation dominated phase. If this is the case, the decay of short wavelength 
modes into superhorizon modes can serve as an active process for seeding superhorizon fluctuations. 

Forthcoming observations of CMB anisotropies as well as large scale surveys with ever greater precision will provide 
a substantial body of high precision observational data which may hint at corrections to the generic and robust 
predictions of inflation. If such is the case these observations will pave the way for a better determination of inflationary 
scenarios. Studying the possible observational consequences of the quantum phenomena found in this article will 
therefore prove a worthwhile endeavor. 
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APPENDIX A: SELF-ENERGY KERNEL FOR k = 



We compute here the kernel /C p (0, rj, r]') 

1 



^p(0; »?,??') = j q 2 dq Im|[5p(g,r?)5|(g,r?')] 2 } , v 



A . 



(Al) 



The infrared divergences at A = | — v = Q arise from the small momenta behavior of the integrand in Eq. l|Al|) . To 
extract such behavior it is convenient to write the Hankel function H^ 1 (qr]) in Eq. I|15|) as [3^| 



[j--„{z)-e l7 < 9 Mz)} 



sax "kv 

The leading q — > behavior follows from the Bessel functions J_ p (z) since [36 



(A2) 



We find upon taking the imaginary part, 

,2 



^lm{[S,(q,r,) S* p ( q> r,')f} q =° 



q 2 



v-X r 2 



r 2 (p) 

2 7T 3 V 



[l + 0{q 2 )} 



(A3) 



The behaviour q 2 ~ 2v = q 2A ~ 1 for q — > in the integrand of Eq. I|A1(I implies a simple pole at A = with a 
77-independent residue (3^|. It is then convenient to add and to subtract from Eq. IjAlfl the low q behaviour Eq. i|A3() . 
We find, 



JCp{0;ri,rf) = h[r],rf',n) + I 2 (r),T]'; 



(A4) 



with the integrals I\ 2 given by 



hM;n) = [" dq (£hn{[S B (q,r,) SUq,v')f} 



q 2 rj rf 



v-X 



2 7r 3 r(i + p) 



' \ v 



2 

rj rj' 



v-X 



r 3 (p) 



2 7T 3 T(l + V) 



dq q 



2-2 V 



1 f°° ( 1 

h{j],r]';n) = — J q 2 dq Im | [S p (q,rj)S*, {q,rj')f | 



(A5) 
(A6) 
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Here, p, > is an arbitrary parameter temporarily introduced to split the integrals. The sum Ii(ry, rj'; fi) + ^2(77, 77'; 7-0 
is clearly /i-independcnt. The second integral in Ii(r),T}'; (j,) [Eq. (|A5() ] can be done straightforwardly, 



„2-2P 



,2 A 



2 A 



(A7) 



This simple integral clearly displays the infrared divergence and the origin of the pole in A: the emission of super- 
horizon quanta for which 577, qq' <C 1 as depicted in fig. [2J 

Keeping the pole in A plus the leading logarithmic terms, and neglecting higher order terms in A, we find for 



1 1 

6tt 2 2A 



rj 2 rj' 2 



+ 



6tt 2 



1 M~H - \ H~m'] 



0(A) 



(A8) 



We can set 3/2 — j/ = A = 0in the integral I2 since it is infrared finite. After lengthy but straightforward algebra, 
we find for prj, fir)' <C 1 



1 

'8^2 



1 



77-77' 



6tt 2 



_n_ _v_ 

7/ 2 r\ 2 



ln[/x(77 - 77')] + 7 + 



1 1 

77' 77 



(A9) 



the principal value P is defined by Eq. Q45J1 • The first term in Eq. \K9i\ is the kernel for a massless conformally coupled 
field (y — 1/2) 27] . It is clear from Eqs. (|A8ll and (|A9|) that the /i dependence cancels out as it should be: 



/Cp(0; 77, 77')=/Ci(0;77,r7') + ^2 
where K, 1 (0; 77, 77') is defined by Eq. 
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2A + 3 
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APPENDIX B: SELF-ENERGY KERNEL FOR fc ^ 



The self-energy kernel /Cp(/c; 77, 77') for the general case is given by Ea. (|21|l . As highlighted in the case of fc = 0, 
there are infrared divergences for v = 3/2 which we regulated with the parameter A = 3/2 — v. We now compute 
/Cp(fc; 77, 77') keeping poles in A and the leading logarithms in 77. The strategy is to separate the regions in the loop 
integral that contain the infrared divergences. Thanks to the azimuthal invariance we can write 



/poo + 1 

d 3 q = 2tt J q 2 dq J d(cos 9) 



with 9 the angle between q and k. Furthermore, we change from the integration variable to p = \q — k\, 



P 



= vV + k 2 - 2 k 



q cost 



d(cos9) = ^- 
k q 



which clearly displays the phase space factor 1/k. Ea. (|21J) takes then the symmetric form 

l-q+k 



Ku(k;r],ri') = — 



q dq 



p dp Im H { ^\qn) H^\qr)') H^iprj) H^iprf) 



r(i) 



q-k\ 



(2)/ 



where we used Eq. (|15H . The integral over p can be performed using [3^ 

J p Hf (pr,) HP (jnf) dp = F B (p, 77, 77') , 



where, 



Fv{p,il,v') = 



rj 2 — rj' 2 



77' H (2) (pr)) Hfliijrq') - 77 H^l^prj) H^iprj') 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 
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We obtain thus for the kernel, 



q dq [F D {q + k, V , r,') - F {\q - k\, V , r,')] H<?\qr,) H^{qr,') 



(B6) 



The infrared divergences at A = I — v = arise from the q — > and q — > k behavior of the integrand in Eq. (|B6|) . To 
extract such divergences we use the small argument behavior of the Hankel functions 38] 



2 7r^ 



<7 2 ?7 77' 



4 
7777' 



[l + 0(9 2 )] , 
g 2 - 2p [l + 0(g 2 ) 



(B7) 



We see from Eas. (|B6ll and (|B7|) that the integrand in Eq.jBGJ) behaves as q = q for q — > and as |g— fc| 



\q — k\ 1 for q — > fc. Therefore, the kernel has a simple pole at A = |39j. It is then convenient to add and to 
subtract from Eq. (|B6|I the behaviour for q — > and (7 — * Eq. (|B7|) as we did in the Appendix A for the k = case. 
We find, 

£p(k;r),r,') = ^ Im jf {<? + M, »/) - - k\,r,,ri')} ttf' \qn) if^W) 



2tt 2 



(F-t:) 



E M 
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(B8) 



v \ 77 77' 



k ImH^'ikrj) Hg } {krf 



(i) 



fi > is an arbitrary parameter temporarily introduced as in Appendix A and we used that 

F,(q + k, v , v ') - F.Qq-kUr,') ^° 2q ^MbM) =2 qk iff (kr,) iifW) • 

/C p (fc; 77, 77') is clearly //-independent as one can easily check by computing the derivative with respect to /1 of the r. 
h. s. of ea. i|B8() . Notice that (9(/i — q) is nonzero for q < /1 while 0(k — q + /i) 9(q — k + fi) does not vanishes for 



k — fi<q<k + (i. The pole at 4 — ^ = A = 0is explicit in the last term of Ea. (|B8|) while the integral over q is finite 
for v = I and fc 7^ 0. 

This analysis for the self-energy kernel shows that all infrared singularites emerge from the regions \qr}\, \qr]'\ -C 1 and 
\prj 1 , \pr)' I <C 1 corresponding to the internal line in the loop which carries momentum q or p = \ q— k | being superhorizon. 
Both regions give a similar contribution because they are equivalent upon re-routing of the loop momentum and Bose 
symmetry. Thus the conclusion of this analysis is that the leading contributions to the self-energy arise from the 
collinear emission of superhorizon quanta [since cos# = 1, see Eq. (|B2Jl j. 

The calculation of the integral for Au(k; 77, 77') in Eq. I|B8|I is straightforward albeit lengthy. These are facilitated 

(1 2) 

by the expression of the Hankel functions H \ ' (qrj) in terms of elementary functions. Dropping contributions of the 

2 

order A the kernel in Eq. I|B8I) becomes for (i<4, \^ir] , \nr]'\ <C 1, 



(B9) 
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(BIO) 



where, 



D(krj,kr)') = —k 3 (w') 3/2 Im 
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=3/2 
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and 
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where, 
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The integral for the kernel takes then form, 
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The remaining calculations are straightforward but tedious. Integrating over q in Ea. l|B12| ) yields for p <C 



/C> /2 (fc; 77,77') 
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47T 2 fc 3 (77 77') 2 
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(B13) 



where a = k(rj — 77') , Ci(z) and Si(z) are the cosine and sine integral functions respectively |38j. The first term in Eq. 
(|B13|) is the kernel for conformally coupled massless particles and the principal part prescription is given by Eq. Q45[l. 

Keeping consistently the leading order terms in A, namely the pole plus finite parts, the dependence on /j, cancels 
out as it should and we find for the kernel to leading order in A 
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(B14) 



Long wavelength limit: |fc?7|, |fc?7'| -C 1 



The behavior of the kernel Eq. (|B13|I in the limit when the wave vector k corresponds to superhorizon wavelengths 
easily follows from Ea. (|B14J) . We use 



Ci[a] Q =° lna + 7 + 0(a 2 ) 



Si[a] a =° a- ^-+0{a 5 ) 
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Gathering all of these results and keeping the lowest order contributions in the limit |A;?7|, \ki]'\ C 1 we find Eq. I|44|) 
for the kernel in the long- wavelength limit. 

Short wavelength limit: |fe77|, \krj'\ 3> 1 
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In this limit the wave function Su(k,rf) is just like the Minkowski space time free field mode function for massless 
fields, namely 



S p (k,J]) = 



^—ikr) 



2k 



(B15) 



It is straightforward to obtain the limit \kr}\, \krf\ 3> 1 of the kernel from Eq. (|B13(I In the short wavelength limit the 
kernel simplifies to 



1 



7] — T} 1 4 7T 2 kr]T\ 

with the constant C given by C = In 2 + 7 — 2. 



7S sinfc(ry- ?/) 



1 ~ , A; 77 77' 
-+C + log^^ 



+ - cosfc(?7- r/) I , 
(B16) 
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